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Abstract 

Some recent work has attempted to show that the singular solutions which are known 
to occur in the Dirac description of spin-1/2 Aharonov-Bohm scattering can be eliminated 
by the inclusion of strongly repulsive potentials inside the flux tube. It is shown here that 
these calculations are generally unreliable since they necessarily require potentials which 
lead to the occurrence of Klein's paradox. To avoid that difficulty the problem is solved 
within the framework of the Galilean spin-1/2 wave equation which is free of that particular 
complication. It is then found that the singular solutions can be eliminated provided that 
the nongauge potential is made energy dependent. The effect of the inclusion of a Coulomb 
potential is also considered with the result being that the range of flux parameter for which 
singular solutions are allowed is only half as great as in the pure Aharonov-Bohm limit. 
Expressions are also obtained for the binding energies which can occur in the combined 
Aharonov-Bohm-Coulomb system. 
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I. Introduction 

In classical physics it is a trivial fact that the absence of a force necessarily implies 
zero scattering. On the other hand the very remarkable Aharonov-Bohm (AB) effect 1 
shows that this does not apply in the realm of quantum mechanics and that potentials 
(as opposed to fields themselves) can indeed have observable consequences. Thus charged 
particles are found to be scattered by a thin magnetized filament even though it is possible, 
by shielding the flux tube or filament, to establish that penetration into the region of 
nonvanishing magnetic field cannot occur. 

For the scattering of a nonrelativistic particle of mass M by the potential 

eAi = a eijTj/r 2 (1) 

where ri is the radius vector in two dimensions and a is the flux parameter one needs to 
solve the Schrodinger equation 

1 (1 



— V — eA J ip = Eip (2) 



2M \i 
Upon writing 

oo 

VM) = jS^/mCr) (3) 



Eq. (2) reduces to the Bessel equation 

Id d 



- - — r -f- + k — (m + a) 2 /r 2 

r dr dr 



fm(r) = ■ (4) 



Since (4) has both a regular and irregular solution, it is necessary to give a boundary 
condition which allows a unique result to be obtained. One could, of course, simply require 
that f m {f) be finite at r = and thereby eliminate ab initio the irregular solution. This 
in fact gives the well known AB solution. Since, however, a resolution of this issue by fiat 
is totally unsuccessful when spin is included, a more physical approach would clearly be 
preferable. This is accomplished 2 by replacing (1) by 3 

eAi = { " C -' ri/r2 T>R (5) 
l0 r < R 
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and taking the limit R — > after matching boundary conditions at r = R. Clearly, the 
vector potential (5) mathematically effects the replacement of an idealized zero thickness 
filament by one of finite radius R which has a surface distribution of magnetic field given 
by 

eH=~8{r-R) . 

Actually, the specific details of the model (5) can be shown to be irrelevant 3 provided only 
that the flux distribution is independent of angle and has no delta function contribution 
at the origin. It is thus straightforward to establish 2 that the irregular solution is absent 
and that the usual AB solution obtains in the R = limit. 

When spin is included, the situation becomes quite different. Here one is concerned 
with the (two-component) Dirac equation 

Eip = [MP + /? 7 • II] V (6) 

where 11^ = — idi — eAi with Ai as in (1). A convenient choice for the matrices in (6) is 

(3 = a 3 

(3(Ji = (ai, sa 2 ) 

where the cr's are the usual Pauli matrices and s = ±1 for spin "up" and spin "down" 
respectively. Upon reduction to a second order form one obtains from (6) the result 



IT 2 + asa 3 - S(r) 
r 



(E 2 - M 2 ) V> = 

which, upon using a partial wave decomposition of the form (3) , becomes 

fm = 



1 d d -2 , , \2 I 2 ! 



- — r — + k — (m + a) Ir — asas - S(r) 
j dr dr r 

where 

k 2 = E 2 -m 2 . 

It is significant that the idealization to a zero radius flux tube has given rise to the com- 
plication that the delta function occurs at a singular point of the differential equation. 
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Different approaches have been attempted in dealing with this difficulty. Alford et 
a/. 4 simply required the upper component to be regular at the origin. This essentially 
makes the spin term trivial and implies an amplitude which coincides with the spinless AB 
result. This contradicts the helicity conservation which is known to be valid 5 for the system 
described by (6). On the other hand, Gerbert 6 has taken an approach which essentially 
states that an arbitrary linear combination of the two solutions J| m + a | (kr) and J_\ m+a \ (kr) 
in the region r > R can be taken so long as it is normalizable for R — > 0. Thus when 
\m + a\ < 1 an arbitrary parameter 9 appears which describes the relative contribution of 
these two functions. This somewhat mathematical approach consequently suffers from the 
appearance of a parameter with no clear physical origin. 

In ref. 3 an approach was taken which was based on the physically reasonable mod- 
ification (5) of the vector potential. Upon matching boundary conditions at r = R and 
letting R — > it was found that the correct solution was always the regular one except 
when as < and |m + a| < 1. In that case there would always be one and only one allowed 
irregular solution. This occurred for 

m = -N , N>0 (7) 

when s = — 1 and 

m = -N-l , N+1<0 (8) 
for s = +1 with the integer defined by 

a = N + /3 

where 

< < 1 . 

It is interesting to note that in the cases (7) and (8) for which irregular solutions are 
allowed there is no contribution from the regular solution J\ m+0 ,\(kr) . Thus the solution 
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obtained in ref. 3 does correspond to a solution of the type obtained by Gerbert without, 
however, the introduction of his mixing parameter. Consequently the solution is unique 
and can be shown 5 to be consistent with helicity conservation. A point worth mentioning 
is that this approach was motivated by a desire to formulate the problem from the outset 
in a physically meaningful way and not (as recently stated 7 ) to provide (after the fact) a 
physical motivation by which to determine Gerbert's 6 parameter. 

It is clear simply from the as < condition that the solution of ref. 3 is not anyonic. 
This has led to some serious concerns by those who consider the anyonic properties of the 
spinless AB system to be of fundamental significance. Thus, for example, ref. 7 has taken 
the introduction of a boundary at r = R one step further by including a nonvanishing 
nongauge potential for r < R and examining its effects when it diverges in certain ways. 
In the following section this work is examined and it is shown that Klein's paradox makes 
that approach unreliable. If, on the other hand, a Galilean spin-1/2 wave equation is used 
in place of the Dirac equation which it closely resembles, then the appearance of Klein's 
paradox can be avoided. This leads to the determination of an "inside" potential which 
allows one to force both regular and irregular solutions to occur. However, it implies 
an energy dependent potential, a fact which is easily seen to follow from dimensional 
considerations. Section III considers the case in which a £/r potential is also included. 
Remarkably, it is found that for arbitrarily small £ the domain of a for which singular 
solutions can be obtained now shrinks to \m + a\ < 1/2. This contrasts sharply with the 
claims of ref. 7, illustrating again the pitfalls of dealing with equations in which Klein's 
paradox is known to occur. 

II. Klein's Paradox and How to Avoid It 

It was proposed in ref. 7 that one modify the spin-1/2 AB problem as done in ref. 3 
by including for r < R a constant potential ur. Since there is no flux for r < R, one has 
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r dr dr 



to consider in that domain the equation 

(E - u R ) 2 - M 2 ] - m 2 /r 2 | / m (r) = 

which has the allowed solution J m (kor) where 

k 2 = (E-u R f-M 2 . 

Since all the results of ref. 7 depend upon the limit of ur — > oo, it is of interest to note the 
behavior of ko as a function of ur. As ur increases from zero ko is a real number (E > m) 
which eventually vanishes at ur = E — m. At this point ko becomes imaginary and the 
function J m (/co^) goes from sinusoidal oscillation to a real exponential (i.e., exp[|/c |r]). 
As ur continues to increase ko vanishes again at ur = E + m and the function J m (kor) 
is again oscillatory and remains so as ur increases without limit. All of these features 
are reminiscent of the phenomenon of Klein's paradox 8 and one concludes that the Dirac 
equation is not an adequate framework for this problem in the limit of arbitrarily large 
Ur. 

On the other hand the type of questions raised in ref. 7 are amenable to treatment 
provided that one substitutes for the Dirac equation its Galilean limit. This is given by 



E 1 - (l + p) + M (l-p)-p 1 .U 



V> = o 



where £ = E — M is the "nonrelativistic" energy. The above was derived in ref. 2 from the 
Dirac equation and is simply the 2 + 1 space version of results obtained by Levy-Leblond 9 . 
It is easily seen to imply for ifti (the upper component of ip) the result 



^i = 



2M 2M R 

provided that one takes a vector potential of the form (5). Again, carrying out a partial 
wave expansion, one obtains the radial equation 



( - ^ + k 2 -[m + a6(r - R)f /r 2 - as)- 8(r - R)\ f m (r) = 

[ r dr dr R J 



(9) 



1 4- — 

1 + \x\ 



has ben introduced for concise- 



where k 2 = 2M£ and the step function 9{x) = \ 
ness. 

One can now include the effect of a repulsive potential as considered in ref. 7. Upon 
letting £ ^ £ — ur one sees that the solution for r < R is J m {kor) with 

/c^ = 2M(£ - . 

At ur = this implies an oscillatory J m (kor) which becomes a real exponential for all ur > 
S. In other words the wave function is exponentially damped in the expected quantum 
mechanical fashion as r decreases from R. The crucial point is that there is no subsequent 
transition back to oscillatory behavior as ur increases without limit and consequently 
Klein's paradox has been eliminated. Thus Eq. (9) is seen to be an appropriate vehicle for 
carrying out the program of ref. 7 which seeks to determine whether a suitably repulsive 
ur can allow the simultaneous occurrence of both J\ m+a \(kr) and J _\ m+a \(kr) for r > R. 
To carry out this study one writes 

f m (r < R) = Im(k r) 

m+a\ 

(kr) + BJ_\ m+a \{kr) 

where I m is the usual Bessel function of imaginary argument and it has been assumed that 
ur > £. The boundary conditions 



Im(\k \R) = AJ\ m+a \{kR) + BJ_\ m+a \{kR) 



dr 



AJ lm+a \(kr) + BJ_ 

\m+a\ 

(kr) - I m (\ko\ 



r 



= asIm(\k \R) 

r=R 



can be solved to yield for A/B the result 

A 



where 





m 


\+g(\k \R) 


2 m + a\ 


+ 


m 


-g(\k \R) 



(fci?)- 2|m+a| (10) 



d 

g{x) = x— log/ TO (x) 



and use has been made of the relation 



\m\ + \m + a\ = —as 

since as shown in ref. 3 only in this partial wave is a singular solution possible. 

One can now answer the question posed in ref. 7 whether a finite nonzero A/B is 
possible. It is, however, clear that since ur — > oo must yield B/A = just as ur = 
gives A/B = 0, there must exist a value for ur which implies a finite A/B. The real issue 
it would seem is whether a ur can be found which is independent of the energy S. It is 
not difficult to show (most trivially, by dimensional considerations) that no such energy 
independent solutions exist. This is in marked contrast to the results claimed in ref. 7. In 
that work the scale for the potential ur is determined by the mass M, or, in other words, 
by a factor which unavoidably requires special relativity. Since Klein's paradox has been 
seen to make that approach unreliable, the scale for ur is determined necessarily by the 
nonrelativistic energy £. 

The (energy dependent) solution of (10) given by 

\k \R = X(kR)\ m+a \ (11) 

where A is arbitrary implies for R — > that 

. /„ 1 A 1 

AB=-- r — . 

4 \m + a\ \m\ + 1 

Thus any (positive) value of this ratio can be obtained by appropriate choice of A. Some- 
what curiously, negative values of A/B can also be generated provided that ur is attractive 
but diverges according to the same power law as in (11). Thus an affirmative answer has 
been found for the issue of fine tuning raised in ref. 7. The fact that this tuning requires 
an intricate dependence of the interior potential on the energy seems to imply, however, 
that it can have little, if any, utility. 
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III. Coulomb Modifications 

It is known that partial wave solutions for the AB problem can be obtained exactly 
even when a 1/r potential is included. In particular such solutions have been obtained by 
Law et al. 10 for the spinless case. This generalization has also been included in ref. 7 in 
the context of their interior repulsive potential ur. Since it has been remarked already 
that Klein's paradox adversely affects such calculations, it is of interest to describe the 
results obtained when the Galilean spin-1/2 equation is employed for such an analysis. 

Upon taking the potential to be 



V(r) 



ur 



r < R 
r > R 



one finds that the appropriate wave equations for the individual partial waves in the 
expansion of the upper component of ip are 



Id d ,o m 

r dr dr r z 



/m(r) = 



for r < R, and 



1 d d , 9 n , rJ . , (m + a) 2 

- ~rr — + k 2 - 2Mi r - ^ 

r dr dr r A 



/m(r) = 



for r > R. It is straightforward to obtain the solution 



fm(r) = J m (k r) 



for r < R, while for r > R 



f m (r) =A rn e lkr (-2ikr^ m+ ^F^\m + a\ + 1/2 + iM£/k 



2\m + a\ + 1 



2ikr 



+ B m e lkr (-2ikrY 



F( - \m + a\ + l/2 + iM£/k 



1 — 2\m + a\ 



— 2ikr 



where F(a\c\z) is the usual confluent hypergeometric function. Note that A m and B m are 
the coefficients of the regular and irregular solutions respectively. It is worth remarking 
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that because of possible modifications which could be induced by the Coulomb term no 
assumptions have been made as to which partial waves can have irregular solutions. 
Upon applying the boundary conditions at r = R there obtains 



^ = L ml (k R)R ^ e ikR (-2ikR)-^F 



\m + a\ + 1/2 + iMijk 



1 — 2|m + a\ 



x-\m+a\ 



F ( - \m + a\ + 1/2 + iM£/k 



1 — 2\m + a\ 



- 2ikrj - e lkH (-2ikR) 
as + R ^jJ\ m \{k Q R) \{e ikR (-2ikRy m+al F( \m + a\ + 1/2 + iM^/k 



2ikR 



2\m + a\ + 1 



- 2ikRj (as + R J\ m \ (k R) - J ]m] (k R)R 
e ikR (-2ikR)\ m+a \F[\m + a\ + 1/2 + iM^/k 



_d_ 

OR 



2\m + a\ + 1 



- 2ikR) 



(12) 



One now takes the R — > limit and finds (as in the ^ = case) that B m must vanish 
unless as < 0, and 

\m\ + \m + a\ = — as . (13) 

However, Eq. (13) is necessary but has not been shown to be sufficient. Nor can one 
merely assume on the basis of the analysis of refs. 3 and 6 that |m + a| < 1. As stressed in 
the former work a second condition emerges when one considers the next-to- leading term 
in powers of R. In the £ = case things are considerably simpler since the solutions for 
both r < R and r > R are Bessel functions whose expansions are characterized by the 
fact that only alternate powers of the argument occur. This is not true for the confluent 
hypergeometric function and one finds from (12) that the expansion in powers of R yields 



B r 



£R 



l-2\m+a\ 



(14) 



for small R. Thus singular solutions are possible for £ ^ only when |m + a| < 1/2 rather 
than the full range \m + a\ < 1 assumed in ref. 7. For the case £ = the result (14) is 
replaced by 

An ^ R 2-2\m+ a \ (15) 

Bm 
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because of the noted property of the Bessel function for small argument. Clearly, the 
condition \m + a\ < 1 follows from (15) whereas the considerably stronger condition \m + 
a\ < 1/2 is required when a Coulomb term is present. It should be remarked that when 
this condition is satisfied, the result (11) again obtains to determine the dependence of ur 
on k and R in the case that both regular and irregular solutions are required to occur. 

The energies S n of the bound states which occur for £ < are readily determined from 
the series expansions of the relevant confluent hypergeometric functions. These yield for 
ur = the results 

i Me 

n " 2 [n-\ ±|m + a|P ' n " ^ ^' ' ' ' 

where the upper and lower signs refer respectively to the case of regular and irregular 
solutions. Of particular interest is the fact that the binding energies of the irregular 
solutions become arbitrarily large as |m + a| approaches \. This illustrates the crucial role 
played by the condition |m + a| < | and lends added credence to the derivation presented 
above. 

IV. Conclusion 

This paper has explored the possibility that more general solutions to the spin- 1/2 AB 
problem can be found. As in the earlier work of ref. 7 this has been done by introducing a 
very short range, but singular, repulsive force. Not surprisingly, it has been found (at least 
in the unambiguous Galilean case) that such solutions can in fact exist. This is physically 
reasonable since the effect of such a potential is to reduce significantly the interaction of 
the magnetic moment with the singular magnetic field at the origin. On the other hand 
such a potential must be required to be energy dependent, and it also has the disadvantage 
of violating helicity conservation, a property which would otherwise be satisfied. Since the 
solution 3 of this problem without the additional non-gauge potential is known to be at 
variance with anyon features, it is clear that those properties can be restored only with 
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considerable difficulty. If one elects to do this, the effect is to negate the full dynamical 
participation of the spin in the interaction. 

One of the most interesting results obtained here has to do with the modifications 
associated with the inclusion of a Coulomb term. It was found that the condition \m + 
a\ < 1 which is generally thought to follow from a condition of normalizability of the 
solution is not sufficient. A more careful analysis shows that only half that range is in fact 
allowed for singular solutions. This has as an immediate consequence that if one considers 
a gas of such particles then the discontinuities known to characterize the second virial 
coefficient 11 B 2 (a,T) are shifted from integer values of a to half-integer values. This is 
particularly noteworthy because of the fact that the transition point has no dependence 
on the strength of the Coulomb potential. Thus one can imagine this parameter to be 
continuously decreased to zero and find that the discontinuities in B 2 (a, T) generally occur 
at half-integers but at integral values when the Coulomb term exactly vanishes. This 
provides a most remarkable example of a system in which a point of discontinuity of a 
variable which has a macroscopic discontinuity has itself a discontinuous dependence on a 
microscopic parameter. This is a subject which clearly merits additional study. 
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